Insulator-to-Metal Transition Induced by Disorder 
in a Model for Manganites 
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The physics of manganites appears to be dominated by phase competition among ferromagnetic metallic 
and charge-ordered antiferromagnetic insulating states. Previous investigations (Burgy et al., Phys. Rev. Lett. 
87, 277202 (2001)) have shown that quenched disorder is important to smear the first-order transition between 
those competing states, and induce nanoscale inhomogeneities that produce the colossal magnetoresistance 
effect. Recent studies (Motome et al. Phys. Rev. Lett. 91, 167204 (2003)) have provided further evidence that 
disorder is important in the manganite context, unveiling an unexpected insulator-to-metal transition triggered by 
disorder in a one-orbital model with cooperative phonons. In this paper, a qualitative explanation for this effect 
is presented. It is argued that the transition occurs for disorder in the form of local random energies. Acting over 
an insulating states made out of a checkerboard arrangement of charge, with "effective" site energies positive 
and negative, this form of disorder can produce lattice sites with an effective energy near zero, favorable for 
the transport of charge. This explanation is based on Monte Carlo simulations and the study of simplified 
toy models, measuring the density-of-states, cluster conductances using the Landauer formalism, and other 
observables. The applicability of these ideas to real manganites is discussed. 

PACS numbers: 75.47.Lx, 75.30.Mb, 75.30.Kz 



I. INTRODUCTION 

The study of manganites - the materials with the colos- 
sal magnetoresistance - continues attracting the attention of 
condensed matter experts.'' 2 These compounds are interest- 
ing for their potential technological applications as "read sen- 
sors" in computers, as well as for the basic science challenges 
that their unusual properties represent for our understanding 
of transition-metal oxides. It is widely believed that strong 
correlations are crucial to understand these manganese ox- 
ides, either in the form of large Coulombic interactions or 
large electron Jahn-Teller phonon couplings, or both simul- 
taneously. Some of their unusual properties include a re- 
markable response to magnetic fields, with the DC resistivity 
changing by 10 or more orders of magnitude in some com- 
pounds upon the application of fields of order 1 Tesla. An- 
other interesting property is the rich phase diagram that these 
materials present, with competition of very different states 
mainly involving ferromagnetic (FM) metallic and charge- 
ordered (CO) antiferromagnetic (AF) phases. Moreover, in 
recent years a plethora of experiments have unveiled a remark- 
able tendency to form inhomogeneities, which occur mainly at 
the nanoscale but with some investigations reporting domains 
as large as submicrometer in size^ These spontaneous gener- 
ation of nanoclusters was predicted theoretically when phase 
separation tendencies were unveiled in the first Monte Carlo 
simulations of models for manganites in 1998 by Yunoki et 
air then confirmed and refined experimentally in dozens of 
effortsji The successful cross-fertilization between theory and 
experiment is remarkable in the area of Mn-oxides, leading 
to considerable progress in recent years. These studies may 
have consequences not only for manganites, but for a vari- 
ety of other compounds where inhomogeneities have been 
found in experiments, including the famous high-temperature 
superconductors m 

Among the leading candidates to understand the CMR ef- 



fect is the 'phase separation' scenarioi 2 ^ In this context, the 
state that is believed to be abnormally susceptible to exter- 
nal magnetic fields is made out of coexisting clusters of two 
competing phases. The clusters involving the FM state have 
random orientations of their magnetic moments, leading to 
an overall vanishing magnetization. However, the "building 
blocks" of the FM state (i.e. the FM nano clusters) are pre- 
formed and for this reason the state reacts rapidly to magnetic 
fields, aligning the nanocluster moments in the presence of 
small fields. Energetically it is quite different to form a FM 
state from "scratch", as opposed to having preformed FM is- 
lands that only need to reorient their magnetic moments to 
generate global ferromagnetism. The competing state is also 
important in the proposed CMR state, since it provides walls 
that prevent the alignment of the nanoclusters in the absence 
of external fields^ 

In order to stabilize this type of CMR states, the role of 
quenched disorder appears to be important. 7 This disorder 
prevents the system from being totally ferromagnetic or totally 
CO-AF, as it occurs in the clean limit where a first-order phase 
transition separates the states^ Small amounts of quenched 
disorder locally favor one phase over the other, leading to the 
inhomogeneous patterns found in simulations, which are be- 
lieved to correspond to those in experiments. These predic- 
tions were beautifully confirmed by the group of Tokura^ by 
means of a careful study of a particular Mn-oxide that spon- 
taneously forms a structure very close to the clean limit, since 
all ions order in a regular pattern. This compounds presents 
a "bicritical" phase diagram, with a first-order FM-AF tran- 
sition in excellent agreement with Monte Carlo simulations. y 
When the material is rapidly quenched in the growing pro- 
cess, such that disorder is now incorporated in the distribution 
of trivalent ions, then a behavior characteristic of many other 
manganites is recovered. In particular, a large CMR is found 
in the presence of disorder, result compatible with the theo- 
retical predictions Recent investigations^ have shown that 
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the addition of lattice elastic effects (due to the cooperative na- 
ture of the Jahn-Teller distortions) makes the disorder strength 
needed to induce the CMR phenomenon smaller in magnitude 
than without elasticity, and avoids limitations in the critical 
dimension that the naive use of Imry-Ma argumentations 11 
would suggest. The relevance of elastic effects has been re- 
marked extensively in recent literature as well Other power- 
law decaying effects, such as unscreened Coulomb interac- 
tions, have also been proposed as relevant in this context^ 
In addition, investigations of a variety of models of percola- 
tion have shown that "correlated" disorder - to mimic elastic 
effects - can even change the order of the transition render- 
ing it first orderpillt may occur that "infinitesimal" disorder is 
enough to trigger and stabilize a phase separation process that 
is intrinsic to systems with first-order phase transitions. Nev- 
ertheless, in practice at least small amounts of disorder appear 
to be crucial to produce the inhomogeneous state that theorists 
proposed as the key factor to understand the CMR physics. 

More recent investigations by Motome, Furukawa, and 
Nagaosa' 5 have provided important information for the the- 
oretical understanding of manganese oxides. Monte Carlo 
simulations of the one-orbital double-exchange model includ- 
ing cooperative phonons have revealed an interesting transi- 
tion from an insulator to a metal upon the introduction of 
quenched disorder into a CO state. The disorder is in the 
form of on-site random energies. Once again, the key role 
of disorder is unveiled by these simulations. Independently, 
Aliaga et aim arrived to analogous conclusions studying a 
two-orbital model for manganites including cooperative Jahn- 
Teller phonons. These authors noticed that disorder drasti- 
cally affects the CE charge-ordered state, transforming it into 
a "CE glass". Overall, it is clear that the use of cooperative 
phonons (for previous work see Refs. Il6ll7l and references 
therein), as opposed to the simpler case of non-cooperative 
lattice distortions, is important to unveil interesting effects of 
manganite models, that may be of relevance to experiments. 

It is the purpose of this work to provide an intuitive under- 
standing of the results of Ref.[3 mainly addressing the tran- 
sition from an insulator to a metal upon the introduction of 
disorder, a result which seems counter-intuitive at first sight. 
Here, it is argued that the transition occurs through an inter- 
esting percolative-like process due to the effect of random on- 
site energies over a charge-ordered state. The sites of the lat- 
tice that had an effective "negative energy" in the CO state - 
namely those that were populated with localized electrons in 
the checkerboard arrangement - may acquire a nearly zero to- 
tal energy due to the influence of that type of disorder. The 
"nearly zero energy" sites can contribute to form a metallic 
state which is highly inhomogeneous in instantaneous Monte 
Carlo snapshots, although it slowly evolves with time (dynam- 
ical) as the Monte Carlo simulation progresses. It is concluded 
that the case found in Refs. I9ll5l provides another very in- 
teresting example where a combination of explicit quenched 
disorder and charge-ordering tendencies leads to nontrivial 
physics, including insulator-metal transitions. 

The organization of the paper is as follows. In Section II, 
the model with cooperative phonons and disorder is discussed. 
The main results of the Monte Carlo analysis are in Section 



III, including the phase diagram, density-of-states with and 
without disorder, and calculation of conductances. A simpli- 
fied model that captures the essence of the main model is in- 
troduced in Section IV. Both models behave very similarly. 
Our proposed explanation of the insulator-metal transition is 
presented in Section V, where the inhomogeneous nature of 
the metallic state is described. Simulations with random hop- 
pings and by freezing the localized spins are also reported in 
this section. Conclusions are presented in Section VI. 



II. MODEL 

The Hamiltonian used in this study contains Mn e g itiner- 
ant electrons coupled to the i2 g Mn core-spins, and also to 
phonons according toili 

H = -* E + h - c -) - 2 j h E Si • Si 

-Xt J~] (Uj-Q - Uj !a )cl a Ci a . (1) 

i,CK,(T 

Here, the first term represents the kinetic energy of the carriers 
hopping between Mn atoms on a <i-dimensional cubic lattice, 
with t the hopping amplitude for the process. denotes 
nearest-neighbor pairs of sites, and the rest of the notation is 
standard. The second term is the Hund interaction between 
itinerant Si and localized Si spins, with Jh the Hund coupling 
constant. The third term describes the energy corresponding 
to the lattice-carrier interaction, with A denoting the electron- 
phonon coupling. Appropriate redefinitions of variables allow 
us to have the energy unit t explicitly in this term (see Reffnj). 
The distortions of the six oxygens surrounding a Mn at site i 
are given by the classical numbers u- h ± a , where in 3D (2D) a 
runs over three (two) directions x, y, and z (x, y). One widely- 
used simplification to the above Hamiltonian, without losing 
essential physics, is to take the limit Jh=oo, since in man- 
ganese oxides the Hund interaction is experimentally known 
to be large. In this limit, the e g -electron spin perfectly aligns 
along the i2g-spin direction. This allows us to reduce the first 
two terms of the double exchange Hamiltonian to;^ 

#DE = - Et T ( Si < S j)4<*J + h - C -]> (2) 
<iJ> 

where 

T(Si, Sj) = -tlcos - cos ^ + sin - sin ?l e ^-*s>]. (3) 
v i) v 2 2 2 2 J 

9\ and <fri are the spherical coordinates of the core spin Sj at 
site i. The new operators d\ create an electron at site i with 
spin parallel to the core spin at i. Within this large Hund- 
coupling approximation, the third term that produces lattice 
distortions becomes: 

#e-ph = -Xt^2(ui- a - U^ a )d\di. (4) 
i.n 
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This tendency is balanced by an opposite crystal elastic energy 
that represents the stiffness of the Mn-O bonds: 



(5) 



In the present investigations, it is important to also consider a 
disorder term of the form: 



(6) 



where Ai is chosen randomly from some distribution (box, 
bimodal, Gaussian) and ni=d\ d\. Hence, the full Hamiltonian 
used in the present study is given by: 



H — i?DE + #e-ph + #ph + #A- 



(7) 



In this paper, we have considered two types of disorder: (i) a 
binary or bimodal distribution where Aj=±A, with A being 
a constant, and (ii) a box distribution where A can take any 
value in some interval [—A', A']. In order to compare the 
results between these two distributions, one should establish 
a relation among them. This can be done by requiring that 
the standard deviation Ax should be the same for both. Since 
(x)=0, one has: 



(Ax) 2 = <x 2 ) 



dx x 2 P(x), 



(8) 



where P(x)=8(x + A) + 5(x — A) for the binary distribution 
and P(x)=l in the interval x G [—A', A'] for the box distribu- 
tion. This gives (x 2 )=2A 2 for the former, and (x 2 )=2A' 3 /3 
for the latter. Hence, the relation between A and A' should be 
A'=(3A 2 )i 

The model used in this paper contains only one orbital in 
the e g sector, instead of the more realistic case of two active 
orbitals. Since using two orbitals increases substantially the 
CPU time, in this first analysis the focus is on the one-orbital 
situation. Also note that Motome and collaborators have stud- 
ied a similar model^ as remarked in the introduction, but they 
implemented the cooperative effect between oxygens (namely 
the effect by which a given oxygen is shared by two octahe- 
dra) in a different manner. This does not represent a problem 
since our results for the phase diagram and density-of-states 
are very similar to those of Refll5l Finally, the technique to 
be used involves the exact diagonalization of the fermionic 
sector, while for the classical degrees of freedom, a standard 
Metropolis is used. Additional details of the numerical meth- 
ods are not presented here since they have been extensively 
discussed in previous literature^ 



III. INSULATOR-TO-METAL TRANSITION INDUCED BY 
QUENCHED DISORDER 

A. Phase Diagram 

The phase diagram of the model discussed in the previous 
section was investigated using spin-spin and density-density 



real-space correlation functions, as well as spin and charge 
structure functions in q-space, which are defined, respectively, 
as follows: 



5 w = ^E< s i- s H 



(9) 



l (x)=^(n i )(n i+x )-l(^(n i )) 2 , (10) 



5(q) = ^(S i -S j )e^( r »- r j), 



n(q)=^(n i )(n j )e^( r »- r j ). 



(ID 



(12) 



In all of the correlations, (O) stands for the thermal average 
of the operator O and is defined by: 



(0) = 



TrIC 



-0H] 



Tr[e-P H ] 



(13) 



where (3 is the inverse temperature, N is the total number of 
sites, and the rest of the notation is standard^ In the follow- 
ing sections, results are given in units of t, unless otherwise 
stated. In addition, all the results shown in this manuscript 
correspond to the case of half-filling n=0.5, where in average 
there is one electron for ever y two sites. This was the elec- 
tronic density studied in Refs. l9ll5l and typically such density 
is easily analyzed numerically (contrary to other cases where 
more complicated charge patterns occur, requiring larger lat- 
tices). Moreover, many real manganites have n=0.5. 

As a criteria to decide whether a particular ordered phase is 
stable, namely whether long-range order is present, the equal- 
time spin and charge correlation functions were used. More 
specifically, a charge-ordered (CO) phase is here defined as 
stable if: (i) the real-space charge correlation n(x) at the 
maximum distance possible in the finite-size clusters studied 
is at least 5% of the maximum value found at distance zero 
n(0). The 5% criteria is certainly arbitrary, nevertheless it is 
qualitatively representative of the physics under investigation 
here. Other criteria, i.e. 10% or 1%, simply lead to phase dia- 
grams slightly shifted up or down in temperature. In addition, 
(ii) maxn(q) must have its maximum value at q=(7r, ir), for 
a checkerboard charge-ordered state to be defined as stable. 
These two criteria give very similar results in practice. For 
the ferromagnetic (FM) phase, a similar definition was used: 
in this case a FM phase is said to exist when S'(x) at the max- 
imum distance is larger than 1% of S(Q), which is already 
normalized to one. In addition, the spin correlation in mo- 
mentum space must have a robust peak at (0,0). The excellent 
agreement found here with previous investigations-'™ leads us 
to believe that our criteria for the determination of the phase 
diagram is robust, and variations around this criteria will not 
change qualitatively the main conclusions. 
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FIG. 1: Monte Carlo phase diagram of the Jn=oo one-orbital model 
as given by the Hamiltonian in Eq.0 for an 8x8 lattice, and overall 
density n=0.5. Results at the values of disorder indicated are shown. 
The error bars are representative of the different critical temperatures 
obtained using different starting configurations, as discussed in the 
text. The number of quenched disorder configurations used for the 
averages varies from point to point, but the smoothness of the results 
indicates that there is no crucial dependence with that number. 



In the actual Monte Carlo simulations, the spin and pho- 
nonic configurations, both assumed to be classical degrees of 
freedom, are updated according to the Metropolis algorithm. 
For the 8x8 lattice, the number of Monte Carlo sweeps per 
site is taken as 10 4 for thermalization, typically followed by 
another 10 4 sweeps or more for measurements. The simula- 
tions were done in a two-fold way. In the first case, the start- 
ing configuration of classical spins and oxygen coordinates 
was chosen to be random, namely a high temperature config- 
uration is selected, followed by a slow cooling down. In the 
second case, simulations started with an ordered stated at low 
temperatures, which is followed by a slow heating up. It was 
observed that both procedures gave fairly consistent results, 
and the fluctuations in the values of critical temperatures are 
contained in the error bars shown. 

Using the criteria outlined in the previous paragraphs, the 
phase diagram of model Eq. Q is shown in Fig. ^ I n the 
clean limit, two phases compete: a charge-ordered phase with 
a checkerboard arrangement of charge (to be expected in the 
large A limit), and a ferromagnetic phase at small A which is 
caused by the double exchange mechanism. 2 The existence of 
these phases and shape of the phase diagram is in excel lent 
agreement with the results of Ref.[3( see a l so Refs. PHI) al- 
though different criteria to define long-range order were used, 
independent programs were written, and even the actual mod- 
els differ in the form in which cooperation is introduced. In 
the clean limit there is a region where both orders coexist at 
low temperature, leading to tetracritical behavior. Other inves- 
tigations of alternative models for manganites produced first- 
order transitions between competing phases with AF and FM 
characteristicsi^2ii£ In the present analysis, an AF state is not 
present and phase competition is probably not so intense as in 
more realistic cases. 

The most interesting results arise upon the introduction of 
disorder in the form of random local energies. In this case, it 
was observed that the FM correlations are not much affected 
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FIG. 2: The criteria used for the determination of the Curie temper- 
ature T c : spin-spin correlation function at maximum distance for the 
cluster studied vs. temperature, for A=0.6 (a) and A=1.2 (b). The 
S(xmax)=0-01 line is 1% of its maximum value, the criteria used 
here to determine FM critical temperatures (indicated by the arrows). 
For Too, n(ir, n) vs. temperature for A=1.0 is shown in (c), for two 
strengths of disorder A, and the approximate critical temperatures of 
the CO-state are indicated. The figures show results for one disorder 
configuration. 



by disorder, but the charge correlations are drastically modi- 
fied. At A in the range 0.6 to 0.8 and low temperatures, a CO 
state in the clean limit becomes a charge-disordered state for 
nonzero A (strength of the on-site disorder). This is sugges- 
tive of an insulator to metal transition. The phase diagram and 
prediction of insulator to metal transition based on the stud 
of the density-of-states was previously discussed in Ref. 
(see also Ref. @) and our investigations confirm their results. 
In addition, studies of conductances discussed in the rest of 
the paper confirm the prediction of an insulator to metal tran- 
sition triggered by disorder. It is interesting to remark that 
the FM state does not seem to play any important role in this 
study, it is the insulating CO-state that is highly susceptible 
to the introduction of disorder. Thus, the actual competi- 
tion at the heart of the phenomenon discussed here appears 
to be among charge-ordered and charge-disordered states, as 
explained more extensively in the discussion below using a 
model where the spin degree-of-freedom is removed. 

Figure [2] illustrates some of the evidence gathered for the 
construction of the phase diagram Fig.^ In parts (a) and (b) 
the spin correlations are shown, at two different values of A, 
as indicated. In both cases, the behavior of the spin correla- 
tions is quite similar before and after disorder is introduced. 
It appears that this degree of freedom is not the crucial for the 
understanding of the phenomena, since it barely changes with 
disorder. More relevant is the study of charge ordering, since 
its associated critical temperature Tco is substantially modi- 
fied at nonzero A. One typical example is shown in Fig. Etc): 
here the charge structure factor at wavevector (%, ir) is shown. 
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FIG. 3: Breakdown of the charge-ordered state with disorder: For 
small temperatures, T=0.10, the CO phase survives the disorder 
[(a),(c) and (e)]. However, at larger temperatures, T=0.30, the initial 
CO state melts with disorder [(b), (d), and (f)]. 



There is a clear difference between the cases with and with- 
out disorder, with a reduction of the critical temperature by 
roughly a factor two. With hindsight based on the discussion 
that is presented in the next section, it is understandable that 
on-site random energies will affect severely a CO state. If 
these random energies are negative at some of the sites that 
were not originally populated in the CO state, likely a trans- 
fer of charge to those sites will be induced. The reciprocal 
occurs with positive random energies in sites that were pop- 
ulated in the clean limit. Thus, the phenomenon presented 
here is expected to be basically related to charge ordering (not 
spin) and also to depend on the particular form of quenched 
disorder used, as confirmed later in the discussion. 

For completeness, results for just one quenched-disorder 
configuration of random energies are shown in Fig. [5] Part 
(a) corresponds to a temperature where both with and without 
the disorder (strength indicated), a CO state is stabilized. The 
charge staggered pattern is robust, even with disorder incor- 
porated. This is caused by the cooperative nature of the lattice 
distortions. On the other hand, at T=0.30 (part (b)) the charge 
staggered pattern is substantially distorted with disorder. This 
fact is also very clearly shown in the charge structure factor 
(part (c)) where a robust (ir, ir) peak is replaced by a fairly 
featureless background at T=0.30 and A=0.4. For each indi- 
vidual configuration of disorder, the charge appears to order 
following a complicated pattern (b). However, the average 
over disorder produces a featureless state. 



B. Density-of-States 

In addition to the equal-time correlations used in the pre- 
vious section to establish the phase diagram, in this effort the 
density-of-states (DOS) was also calculated. The procedure to 
calculate dynamical quantities involves the full diagonaliza- 
tion of the electronic problem for each configuration of clas- 
sical spins and phonons. Details can be found in Ref.yJ 

Interesting results are shown in Fig.|4] illustrating one of the 
main conclusions of the present investigation, and confirm- 
ing the predictions of Ref. [15J As expected, the clean-limit 
case A=0 presents a DOS which has a gap at the chemical 
potential, in agreement with the results of Fig. ^ that locate 
A=0.7 and low temperature within the CO phase. The pres- 
ence of oscillations at nonzero values of u-fj, is due to well- 
known finite-size effects: a relatively small cluster with pe- 
riodic boundary conditions typically has a shell structure in 
their ladder of states. These spurious oscillations do not af- 
fect the conclusions of our study. The main result related with 
Fig.|4]is the interesting dependence of the DOS with the mag- 
nitude of A. In fact, at A=0.3 and higher the original CO 
gap at the chemical potential has simply vanished, suggesting 
the instability of the insulating CO-state toward a novel state 
(presumably metallic) induced by quenched disorder. Under- 
standing the origin of this unexpected insulator-to-metal tran- 
sition, triggered by disorder, will occupy most of the rest of 
the present paper. 
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FIG. 4: Density-of-states obtained with Monte Carlo simulations at 
the values of A and temperature indicated. The lattice is 8x8 in size, 
and averages over five disorder configurations are presented (the dis- 
order was generated with a bimodal distribution in this example). 10 4 
thermalization steps were used, with 10 4 measurements. Values of A 
are indicated. Similar results were also obtained using 4x4 clusters 
(not shown). 

To convince the reader that the apparently small number 
of disorder realization used in Fig.@]is sufficient for our pur- 
poses, in Fig.[5]results for five independent quenched-disorder 
configurations are shown. Qualitatively, the five present states 
at the chemical potential as observed in the average Fig.0] 

To analyze the dependence of the results with the type of 
disorder introduced, in Fig.[6]results for a box distribution, in- 
stead of bimodal, are shown at the couplings indicated. Once 
again, the DOS presents a CO state at small A, which is re- 
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FIG. 5: First five figures correspond to the DOS obtained with five 
different quenched-disorder configurations, at the same A, lattice 
size, and temperature used in Fig.[4] A bimodal distribution of disor- 
der was employed and results for A=0.4 are shown. All five results 
are qualitatively similar. The last figure is the average over all five. 



placed by a metallic-like DOS when A is increased. 

To study the dependence of the results with A over a wider 
range, in Fig.Qresults for A as large as 1.6 are shown, us- 
ing a bimodal distribution. At large A, it is natural to expect 
that the states at the Fermi level will not be filled, since this 
strong disorder must induce localization of charge at the sites 
where the on-site random energy happens to be large and neg- 
ative. Then, the DOS presents an interesting non-monotonic 
behavior with an insulator at both A small and large, while 
the system is metallic at intermediate A. 
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FIG. 6: Density-of-states (averaged over five disorder configura- 
tions) using a box distribution (instead of bimodal). 10 4 thermal - 
ization steps were used, with 10 4 sweeps for measurements. 

The situation is different for the box distribution. In this 
case, it is expected that large A will not alter the results so 
dramatically as for the case of the bimodal distribution, since 
there are always random energies that will happen to be small 
and they do not lead to localization. This is confirmed by 
the results of the simulation in Fig.|S] showing that states are 
present at the Fermi level even for large disorder. 



C. Landauer Conductance 

While the presence of weight at the Fermi level in the DOS 
is suggestive of a metallic state, it may occur that those states 




FIG. 7: The "closing" and "re-opening" of the gap in the density- 
of-states as the disorder strength A is increased, using a bimodal 
distribution. Temperature, A, and A's are indicated. 




FIG. 8: Analogue of Fig.|7|but for a box distribution. In this case, 
the gap that closes at intermediate A does not form again for larger 
values of the disorder strength. 



are not extended but localized. To clarify this issue the con- 
ductance, G, of clusters described by the Hamiltonians used 
in this effort was also evaluated. The calculation was carried 
out using the Kubo formula adapted to geometries usually em- 
ployed in the context of mesoscopic systems^ The actual ex- 
pression for G is 



G 



e 2 
2— Tr 
h 



(ihv x )lTaQ(Ep)(ihv x )IiaQ(Ep) 



(14) 



where v x is the velocity operator in the x direction (assum- 
ing current flows along that direction) and ImQ(E F ) is ob- 
tained from the advanced and retarded Green functions using 
2ilmg{E F )=g R (E F ) - G A (E F ), where E F is the Fermi en- 
ergy. The cluster is considered to be connected by ideal con- 
tacts to two semi-infinite ideal leads, as represented in Fig. [9] 
The current is induced by an infinitesimal voltage drop. This 
formalism avoids some of the problems associated with fi- 
nite systems with periodic boundary conditions, such as the 
fact that the optical conductivity is given by a sum of Dirac 
S functions due to the discrete nature of the matrix eigenval- 
ues. Finding a zero-frequency Drude peak with finite weight 
corresponds, in principle, to an ideal metal - i.e. a system 
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with zero resistance - unless an arbitrary width is given to 
the zero frequency ^-function. In addition, in numerical stud- 
ies sometimes it occurs that the weights of the zero-frequency 
delta-peak are negative due to subtle finite-size effects^ For 
these reasons, calculations of DC resistivity using finite closed 
systems are rare in the literature. All these problems are 
avoided with the formulation described here, which can be 
readily applied to the Hamiltonians employed in our study 
where the fermionic sector is quadratic. The problem basi- 
cally amounts to the calculation of transmissions across non- 
trivial backgrounds of classical spins and lattice distortions. 



left 
lead 




right 
lead 



e-oA=0.0 
q-qA=0.4 




=S3i 



0.08 0.12 
Temperature (t) 



0.20 



FIG. 11: Logarithm (in base 10) of resistivity (inverse of conduc- 
tivity) vs. temperature for the same lattice and couplings used in 
Fig. 1101 Note the similarity of the results with those obtained vary- 
ing magnetic fields near AF-FM first-order transitions (see Ref.[<§). 



FIG. 9: Geometry used for the calculation of the conductance. The 
interacting region (cluster) is connected by ideal contacts to semi- 
infinite ideal leads. The configurations of classical localized spins 
and phonons are produced by Monte Carlo simulations. 

In practice, the entire equilibrated cluster, as obtained from 
the MC simulation, is introduced in the geometry of Fig. [9] 
The ideal leads enter the formalism through self-energies at 
the left/right boundaries, as described in Ref. In some 
cases a variant of the method, explained also in Ref. 
was used to calculate the conductance: Instead of connect- 
ing the cluster to an ideal lead with equal hoppings, a replica 
of the cluster was used at the sides. This method, although 
slightly more CPU time consuming, takes into account all of 
the Monte Carlo data for the cluster, including the periodic 
boundary conditions. The final step in the calculation is to 
carry out averages over the localized spin configurations and 
oxygen coordinates provided by the Monte Carlo procedure. 
The physical units of the conductance G in the numerical sim- 
ulations for 2D are e 2 /h as can be inferred from Eq.[2] In 2D 
the conductivity a coincides with the conductance, while in 
3D a is simply given by the conductance divided by the linear 
size of the lattice, assumed cubic. 
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FIG. 10: Conductance (e /h) vs. temperature for the one-orbital 
model with cooperative phonons at A=0.7 and using a 12 x 12 cluster, 
in the geometry shown in Fig.|9l 



A representative result of our effort is shown in Fig. El 
where G vs. temperature is shown at an electron-phonon cou- 
pling and disorder strength where the DOS study of the pre- 
vious section revealed the insulator to metal transition. In the 
clean-limit, the nearly vanishing value of G at low tempera- 
tures reveals the presence of an insulator, which it is known to 
be a charge-ordered state. The transition to a disordered state 
occurs with increasing temperature at approximately T=0.08- 
0.10 (A=0), value compatible but slightly lower that found 
through the charge correlations (this may simply be due to 
the different lattice sizes used). The most remarkable result 
is shown at A=0.4. Here, the conductance is very similar to 
the clean-limit result at temperatures above the transition, but 
at low temperatures they are drastically different. Compati- 
ble with the finite DOS at the Fermi level found in the previ- 
ous section, the system with nonzero quenched disorder has 
a finite conductance at low-T, which actually increases with 
decreasing temperature as in a metallic state. The inverse of 
the conductance (i.e. the resistivity in 2D) vs. T is shown in 
Fig. using a logarithmic scale to amplify the nearly van- 
ishing conductance of the insulating state of Fig. El The in- 
troduction of quenched disorder transforms an insulator in the 
clean limit to a metallic state in the dirty case. As a conse- 
quence, the states at the Fermi level induced by disorder that 
were observed in the previous section indeed correspond to 
extended states that transport charge, a very interesting result. 
Note, however, that the conductance in the metallic state is 
substantially smaller than its maximum value (10 for a 12x12 
cluster), indicative of a dirty metal. This result is also compat- 
ible with experiments in manganites that have revealed simi- 
lar characteristics. 1 - 2 Then, the metallic state is far from per- 
fect and there must exist sources of scattering that prevents 
the smooth flow of charge. In subsequent sections, it will be 
argued that this metallic state is percolative-like, result com- 
patible with its relatively high resistance. 

In Fig. results for the conductance at low temperature 
as a function of A are shown for two different cluster sizes. 
The size effects appear to be very small, although the clusters 
studied are not sufficiently large to show convincingly that a 
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FIG. 12: Conductivity vs. disorder strength for model Eq.0 at A=0.7 
and T=0.05. 10, 000 thermalization and 2000 measurement Monte 
Carlo sweeps were used. Results for two lattice sizes are shown, 
suggesting that finite-size effects could be small. 



metallic state was indeed stabilized. The behavior as a func- 
tion of A reveals a range of insulating behavior at small A, 
followed by a well-defined peak in the conductance at inter- 
mediate values - where A is comparable to the gap, as ex- 
plained later in the text - followed by a subsequent decrease 
of G at A of order 1 or larger. This large- A behavior can eas- 
ily be understood due to strong localization of carriers in that 
limit. However, the presence of a metallic state at interme- 
diate A - as opposed to a smooth interpolation between the 
two insulating limits at small and large A - is puzzling. Its 
understanding is the main goal of the next section. 



IV. RESULTS FOR A SIMPLIFIED MODEL 

The results presented in the previous section have revealed 
an insulator to metal transition induced by on-site random- 
energy disorder, when the clean limit starting state is charge 
ordered. This conclusion appears clearly both in the DOS and 
in the conductance, and it is already evident in lattices as small 
as 4x4. Then, the mechanism must be easy to understand 
since it cannot depend on subtleties related with the bulk limit 
or large correlation lengths. In order to clarify the origin of the 
insulator-metal transition, here a "toy model" involving only 
electrons will be introduced, and shown to behave quite sim- 
ilarly as the actual "real model" that has electrons, localized 
spins, as well as lattice degrees of freedom. The definition of 
the model is: 



(U) 



h.c.) + ^ain ; A i 



(15) 



This simplified toy model includes a hopping term for spin- 
less fermions (first term), as in the realistic case at large Hund 
coupling. The hopping amplitude iy will first be considered 
constant and equal to t in this section, but later a random hop- 
ping will also be studied for completeness. The most impor- 
tant simplification of Eq. [2] as compared to Eq. Q is the 
following: instead of lattice distortions to induce charge or- 
dering, a mere modulation of the chemical potential is intro- 



duced "by hand" (second term in Eq.I15>. where ai=+a and -a 
are located on a checkerboard arrangement. This is sufficient 
to generate the analog of the CO state induced by phonons, 
namely a charge modulation leading to a checkerboard pat- 
tern. In all the toy model simulations reported in this paper 
a was kept equal to 0.5: for computational simplicity it was 
not modified to mimic different values of A and concomitant 
charge gaps. However, the qualitative results emerging from 
the analysis are sufficiently clear, that there was no need to 
further consider a as a tuning variable. Finally, to incorporate 
the quenched disorder, an on-site random-energy term as in 
the realistic case was included as well. 



1.00 




1.8 2.0 



FIG. 13: Conductivity vs. disorder strength A, comparing results for 
the "real" phononic model Eq.0at A=0.7 against the "toy" model. In 
both cases a bimodal distribution is used, and results are shown using 
a 12x12 lattice at T=0.05. In the calculation of the conductance 
G for the real model, 10, 000 sweeps for thermalization and 2000 
for measurements were used. For the toy model the curve shows 
an average over 100 configurations of disorder, and a is 0.5. The 
different locations of the conductance peaks can be brought to the 
same value by modifying a, an unnecessary task in our qualitative 
study. 

The calculation of conductances vs. disorder strength, as 
well as other observables, proceeds straightforwardly for the 
toy model as a special case of the realistic one. There is 
no need to carry out any Monte Carlo simulation for the toy 
model, the results are exact. A remarkable result is shown in 
Fig. where G vs. A is shown for both the realistic and 
the toy models, on the same lattice and at the same temper- 
ature. Besides a harmless shift in the location of the peak 
feature, the shape of the curve is basically the same, revealing 
insulator-metal-insulator characteristics with increasing A in 
both cases. Then, the simplified toy model is sufficiently ro- 
bust to have the same properties as the real model, and its 
analysis will likely unveil the origin of the effect. This is con- 
firmed by a study of the DOS as well (Fig. ll4b . since the gap 
in the clean-limit caused by the staggered-modulation in the 
local-energy induced by the second term of Eq.^]is filled by 
the introduction of disorder. 

The analysis of the real model - i.e. including phonons and 
localized spins as active degrees of freedom - was limited to 
12x12 lattices and a few realizations of disorder. With the toy 
model, a far better study can be conducted. For example, on 
Fig-El an analysis of G vs. A at several lattices is shown. 
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FIG. 14: Density-of-states of the toy model Eq. ll5l for several values 
of the disorder strength A. The lattice is 16 x 16 and temperature 
T=0.05. At the Fermi level states are generated, as it occurs in the 
realistic model studied in previous sections. 



All of the curves have the same shape, and the size effects are 
moderate. However, there is a slow decrease in the conduc- 
tance at the peak. Given the error bars of the simulation (the 
oscillations in the results between subsequent A's can provide 
a crude estimation of these error bars), and the rapid growth 
of CPU time with increasing size, an even more careful analy- 
sis cannot be carried out. As a consequence, it is not possible 
in this early investigation to be conclusive about the metallic 
nature of the state in the bulk 2D limit, particularly in view of 
the subtleties that are sometimes related with similar studies 
in the case of Anderson weak-localization effects (this partic- 
ular issue is quantitatively discussed in the Appendix). Nev- 
ertheless, the qualitative analysis to be discussed below will 
be sufficient to understand the nature of the metallic state in 
finite-size clusters, as a first step toward applications to real 
compounds. 
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FIG. 15: Conductivity vs. disorder strength for the toy model, aver- 
aged over 100 disorder configurations at T=0.05. Results for several 
cluster sizes are shown. In 2D, conductivity and conductance are the 
same. The size effects appear to be mild, but the slight reduction of 
G with lattice sizes in the central region could still lead to an insulat- 
ing state in the bulk limit. More work is needed to fully understand 
this subtle issue. 

The situation in 3D is similar, since the results for the con- 
ductivity (G/L, where the lattice size is L 3 ) of the toy model 



seem to converge with increasing size of the cube used for the 
simulation, as shown in Fig. [H)] While it is not impossible 
that the conductivity could slowly decrease to zero as L^oo, 
the results are suggestive that a metallic state could indeed 
be found in the bulk 3D limit. More work is needed to fully 
confirm this assumption. 
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FIG. 16: Conductivity (in e 2 /hL units) vs. disorder strength for the 
toy model in 3D, averaged over 100 disorder configurations. L is the 
size length of a L 3 cube. 



V. UNDERSTANDING THE METAL-INSULATOR 
TRANSITION 

A. Proposed Explanation 

The essence to understand the insulator to metal transition 
found with increasing disorder-strength appears to be very 
simple, and a starting point is illustrated in Fig.^] Suppose 
that in the clean-limit case, A=0, the Fermi level is in between 
two states separated by a gap (and, as a consequence, the state 
is an insulator). It is natural to expect that with increasing on- 
site energy disorder those two states will increase their respec- 
tive widths - since the on-site random energies directly affect 
the energy location of the state - and there will be a value of 
A (of the order of the clean-limit gap) where a robust number 
of states will be located at the Fermi energy. The results we 
have found in our investigation are compatible with this sim- 
ple picture for the generation of weight at the Fermi level. In 
this respect, increasing temperature in the clean-limit at fixed 
large A or increasing A at fixed small T and large A appear 
to be analogous procedures to transform an insulator into a 
metal. However, Fig. [^]can only be a rough starting point 
to rationalize the metallic conductance since the presence of 
states at the Fermi energy is not sufficient to guarantee that the 
systems is indeed metallic. In fact, Fig.[n]applies even in the 
limit of vanishing hopping, where the charge is localized. 

The most revealing procedure to understand the I-M-I tran- 
sition with increasing A found in the real and toy models in- 
volves the study of Monte Carlo "snapshots" of the charge 
configurations. This is useful in the real model because 
the classical degrees of freedom freeze at low temperatures, 
and snapshots are representative of the physics (a state with 
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FIG. 17: Cartoonish view of the generation of weight at the Fermi 
level by disorder. If in the clean limit the system is insulating (with 
no states available at the Fermi level), the mere broadening of these 
states by disorder certainly introduces states at Ef after some A 
strength similar to the gap is reached. However, it remains to be in- 
vestigated whether the states at the Fermi level generated by disorder 
are truly metallic, or whether they are insulating. 

1=0.7 l=0J h=0.7 
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FIG. 18: Monte Carlo "snapshots" obtained using a 12x12 cluster 
and temperature T=0.05 showing the electronic density at each site 
in tones of black and white. Results are shown for both the "real" 
and "toy" models (Eqs.l7land ll5l respectively). In the former, A=0.7 
is used. In the three cases, results for the clean limit A=0, as well 
as intermediate and large values of A are presented. Both the real 
and toy models at A=0 have densities approximately 0.7 and 0.3 in 
a checkerboard pattern, while at A=2 the two clearly-distinguished 
densities are close to 0.95 and 0.05. A detailed discussion is in the 
text, but the reader can already visualize the localized nature of the 
charge in both limits of small and large A, while the snapshots at 
intermediate disorder show a more disordered distribution of popu- 
lated sites, leading to a metallic state. It is also to be remarked the 
similarity between the results for both models. 



strong quantum fluctuations may have many configurations 
with equal weight, complicating the analysis). Typical snap- 
shots are shown in Fig. where results for the electronic 
density of both the real and toy models are contrasted. It is 
clear that in the three regimes of relevance - small, intermedi- 
ate, and large A - there is an excellent agreement between the 
two cases, and the toy model appears to capture the essence 
of the behavior of the more realistic model Eq.0 Let us now 
analyze the results. At small A, the checkerboard pattern cer- 
tainly prevents the transmission of charge, and the system is 
insulating in both cases. At large A, the sharp color-contrast 
between the sites occupied and those empty are indicative of a 



strong localization of charge. Energetically this occurs to take 
advantage of the regions with large and negative A, where at 
every site either occupancy zero or one is favored. This regime 
is insulating as well. 

The most interesting result occurs at intermediate values of 
A, in the metallic regime. Here the density is closer to its 
mean value 0.5 in several lattice sites, in spite of the a modu- 
lation of the toy model. To understand the results, consider for 
instance the case of a bimodal distribution for the quenched 
disorder, as in Fig. ED Considering the second and third terms 
of Eq.^]in the toy model, then locally the effective chemi- 
cal potential can take four values, namely A+a, A-a, -A+a, 
or -A-a. The case of interest for metallicity is obtained at 
A~a or larger - in excellent agreement with the numerical 
results - where approximately half the sites have a small ef- 
fective on-site energy, and these sites are now available for the 
transport of charge with average density 0.5 (while the other 
sites either have a very large or very small chemical potential, 
favoring either minimum (zero) or maximum (one) local oc- 
cupancy, with both cases not useful for charge transport). The 
sites with small effective on-site energy would not contribute 
to a metal if they were isolated (for instance if they were reg- 
ularly spaced), but for a random distribution of locations a 
metallic state can be formed if there is a large enough num- 
ber to induce percolative-like phenomena. Figure [T9l further 
clarifies this issue. The calculations of conductances appear 
to imply that indeed such percolation has occurred at the den- 
sities here considered. This appears to be the essence of the 
insulator to metal transition. The state that conducts appears 
to be highly inhomogeneous, compatible with the low value 
of its conductance. Clearly, insulating regions survive. The 
elastic cooperative nature of the oxygen coordinates enlarges 
the regions that are either charge ordered or disordered. 



(a) (b) 




FIG. 19: Example showing that the sites with densities approxi- 
mately 0.5 are correlated with those with small effective on-site en- 
ergy, for the particular case A=0.8, and a=0.5, on a 24x24 lattice. 
On the right panel, sites with small effective on-site energy are shown 
in black if A-a=0.3 and light grey if A-a=-0.3, the rest are in white. 
On the left panel, results for the toy model are shown. Here sites with 
densities between 0.50 and 0.75 are in black, between 0.25 and 0.50 
in light grey, and the rest in white. The agreement between left and 
right is almost perfect. This two-colors view helps in the understand- 
ing of percolation. In two dimensions the critical percolative fraction 
is 0.5, precisely the number of sites with small effective on-site en- 
ergy. The system is then only at the verge of percolation. However, 
in 3D it would be fully percolated since the critical fraction is 0.25. 
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To clarify further the poor nature of the metallic state, a 
frequently-mentioned issue in the study of conductances is 
their distribution. Sometimes the mean value, as the one dis- 
cussed thus far, is not representative due to the large width of 
the distribution. However, an example of G-distributions is 
shown in Fig.[20| for the toy model at a value of A where G 
indicates a metallic state. The result shows that the distribu- 
tion is not abnormally wide, but the only values that visibly 
appear are 0, 1, and 2, while its full possible range of values 
extends up to 22 in the case of the largest lattice studied. The 
distribution shown also indicates that slightly less than half the 
disorder configurations used have zero conductance, while the 
others have the minimum and next values, borderline compat- 
ible with a metal. Then, the metallic state discussed in this 
paper corresponds to a poor metal, and the snapshots previ- 
ously discussed indicate that this property arises from inho- 
mogeneities. The system appears to be near percolation and, 
as a consequence, there are only a limited number of paths for 
conductance across-the-sample. 




FIG. 21: Density-of-states for the toy model where now the disorder 
is introduced through random hoppings, rather than random on-site 
energies. Here, the hopping disorder is again a bimodal function 
with values ±r. The results are at T=0.05 and obtained on a 12x 12 
cluster. The gap does not close under the influence of disorder in the 
hoppings. 
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FIG. 20: Histogram of conductances observed in the study of the 
toy model, corresponding to A=0.9, and at T=0.05, for many values 
of the disorder. The conductances can only take integer values, here 
they have been broadened for a better illustration. Size effects appear 
to be small. Several configurations are insulating (G=0), while others 
conducting (although poorly). 



Although the use of standard formulas for percolative pro- 
cesses could be tempting (for instance to predict the critical 
density of metallic sites needed to percolate, and the corre- 
sponding critical exponents), here the electronic-density at 
each site is not simply found to correspond to just four num- 
bers, as the effective on-site energies argumentation would in- 
dicate. The reason is that the kinetic energy of electrons (first 
term in Eq. 1 151 must also be considered, and the electronic 
density "spills" from favorable to unfavorable sites if they are 
close to one another. This proximity or spilling effect is not 
considered in traditional models of percolation - where links 
or sites are sharply either metallic or insulating - preventing 
the straightforward use of standard formulas deduced for clas- 
sical random-resistor networks to the case analyzed here. 



B. Consequences of the Proposed Explanation 



1. Results with Random Hoppings 



The model used both here and in the previous investiga- 
tion^ employs a local random-energy term as a source of dis- 
order. This form of disorder appears to be crucial for the argu- 
mentation discussed in the previous section to explain the sta- 
bilization of a poor metallic state at intermediate A. The no- 
tion of an effective site-energy that nearly cancels at A^a for 
approximately half the sites is very important to render those 
sites conducting. Then, the present analysis would predict 
a lack of universality of the disorder-induced insulator-metal 
transition upon modifications in the explicit form in which the 
disorder is introduced. In particular, introducing randomness 
in the hopping amplitudes, rather than in the on-site energies, 
should not lead to a metal according to the previously pro- 
posed explanation of the effect, since the effective site ener- 
gies will not be affected. To verify this idea, in the toy model 
the on-site disorder was switched off, but now the hopping 
amplitudes were given an extra random component. More 
specifically, the hoppings were modified such that ijj=i+Tij, 
while Aj was switched off in Eq. ^] The random number ry 
was chosen from a bimodal distribution with values ±r. 

The results for the density-of-states are shown in Fig.l21lfor 
several values of r. There is a drastic difference when com- 
pared with the case of random on-site energies. For random 
hopping amplitudes the gap present at r=0 is not filled by the 
random hoppings with increasing r. The state remains insu- 
lating at all values of t, small, intermediate, and large. This is 
compatible with the explanation for the insulator-metal tran- 
sition discussed in the previous subsection, which relies cru- 
cially on the effect of on-site disorder. 
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FIG. 22: Phase diagram of model Eq.0obtained in the limit where 
the localized spins are frozen in a ferromagnetic state at all temper- 
atures. As in the case of non-frozen spins, 10 4 thermalization and 
10 4 measurement Monte Carlo sweeps were typically used at each 
temperature and A. The lattice was 8x8. The criteria used to deter- 
mine the transition temperatures are the same as for the simulations 
with active spins. Results both with and without disorder are shown, 
with a behavior similar to that found in Fig.Q The phase diagram at 
very low temperatures is difficult to obtain due to the rapid increase 
in CPU time required to remain ergodic (configurations are separated 
by large energy barriers). 



2. Simulations with Frozen Spins 

Another important prediction of the proposed explanation 
for the I-M transition is the irrelevancy of the spin degree-of- 
freedom, and its concomitant FM phase. If the idea is cor- 
rect, then the presence of the CO vs. FM phase competition 
is not crucial, but more important is the competition between 
charge-ordered and charge-disordered states. To verify this 
idea a simulation was carried out freezing the localized spins 
to a perfectly ferromagnetic state at all temperatures, effec- 
tively removing the spin as a degree-of-freedom. The Monte 
Carlo phase diagram is shown in Fig.[23]and indeed it shows 
the expected behavior for the CO state, before and after in- 
cluding disorder. 

The density-of-states for the model with only the charge as 
active degree-of-freedom is shown in Fig. [23] with and with- 
out disorder. The insulator at A=0 and metal at A=0.4 are 
present qualitatively as in the case with active spins. It ap- 
pears that the insulator-metal transition only depends on the 
presence of charge ordered and disordered states, and their 
competition. The ferromagnetic component is of no conse- 
quence for the effect discussed in this paper. 



VI. CONCLUSIONS 

In this paper, an explanation was proposed for a recently 
discovered insulator-to-metal transition induced by disorder- 
in a model for manganites with cooperative phonons. The 
transition occurs when the disorder is introduced in the form 
of random on-site energies. Our effort presented here confirms 
the effect, showing that the clean-limit gap in the density-of- 



FIG. 23: Density-of-states obtained using the model with frozen 
spins, leaving only the charge as active degree-of-freedom. The 
clean-limit gap is still found to close with disorder, even in the ab- 
sence of CO-FM competition. The lattice is 8 x 8. 



states is filled by disorder. Moreover, a calculation of the 
conductance was implemented, with a Landauer formulation 
borrowed from the context of mesoscopic physics. This cal- 
culation confirms that for intermediate values of the disorder 
strength A the system is indeed metallic, at least for the clus- 
ters studied here. A better finite-size scaling should be carried 
out to fully confirm the metallic nature of the state in the bulk 
limit, but this huge effort is left for future investigations. Our 
limited size analysis thus far indicates that the conductivity 
will remain finite in the bulk limit, particularly in three di- 
mensions. 

To understand the origin of this transition, an even simpler 
model without active phononic degrees-of-freedom was pro- 
posed. This new model behaves very similarly as the original 
one, with a DOS that is filled by disorder, a quite similar con- 
ductance vs. A curve, and similar Monte Carlo snapshots. 
Both this model and the original one, present a percolative- 
like behavior in the metallic region at intermediate A's. The 
reason is that there are sites where the random energies com- 
pensate a dynamically generated on-site energy induced by 
the interaction with phonons. This interaction produces a 
checkerboard pattern of charge, leading to sites that are either 
charge localized n=l (negative effective energy) or empty n=0 
(positive effective energy). The compensating random energy 
induces some of the sites to prefer charge n~0.5, as in a ho- 
mogeneous metal. It is through these sites that charge can 
move from one side of the cluster to the other. Compatible 
with this picture - where only a fraction of the sites have n 
approximately 0.5 - the metallic state is found to have poor- 
metallic behavior, and in the histogram of conductances for 
different realizations of quenched disorder many of them have 
insulating character. The mechanism here discussed to un- 
derstand the filling of gaps by disorder appears general, and 
it should be valid beyond manganites. In fact, there are ex- 
perimental results and theoretical investigations that have re- 
ported similar effects in the context of disordered insulators 
(see Ref.|^l and references therein). In these related investi- 
gations random configurations of disorder were found to pro- 
vide conducting resonant electronic channels for tunneling in 
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tunnel junctions, an explanation similar to that provided in the 
present investigations. In addition, it is known that several ma- 
terials share phenomenological aspects with the manganites, 22 
and the investigations of these interesting effects may have 
consequences in a variety of research areas. 

Is this mechanism active in the real manganites? More work 
is needed to investigate this issue. Reasons for concern are the 
dependence of the results with the form of disorder, and the 
irrelevance of the competing FM state in the process. Disor- 
der in Mn-oxides is believed to arise from lattice distortions 
induced by chemical doping with ions that are either much 
larger or smaller than those involved in the parent undoped 
compound. From this perspective disorder in the hopping 
amplitudes would be realistic. On the other hand, the ionic 
dopants have often a different valence than the ion they re- 
placed and that could be represented by on-site energies. Ad- 
ditional effort is needed to relate this simple insulator-metal 
mechanism to those proposed in early investigations and to 
experiments. However, even if the immediate relevance of the 
idea to Mn-oxides needs further study, the effort of Refs. l9ll3 
and ours confirms the key role that disorder plays to render 
metallic an insulating system, a somewhat counter-intuitive 
idea. This adds further evidence that the most basic aspects of 
the original scenario by Burgy et ali are essentially correct, 
namely disorder is needed at least as a triggering effect to in- 
duce insulator to metal transitions in models for manganites 
in the vicinity of regions where FM-AF phase competition oc- 
curs. Moreover, the present effort has shown that the metallic 
state induced by disorder is "poor" due to its inhomogeneous 
characteristics. All these results provide further support for 
the mixed-phase-based percolative explanation of the remark- 
able DC transport behavior of the CMR manganites. 
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VIII. APPENDIX: ISSUES RELATED WITH ANDERSON 
LOCALIZATION 

It is clear in our figures for the conductivity in 2D and 3D, 
corresponding to the simplified "toy" model, that a (or the 
conductance) reaches its maximum value at some intermedi- 
ate disorder strength A^O.8-0.9 (see Figs. I15landll6i. How- 
ever, both in 2D and 3D the conductance maximum value de- 
creases slowly as the lattice size increases. As discussed in 



the text, this prevents a clear statement regarding the stability 
of a metallic phase in the bulk limit by the mechanism dis- 
cussed in this paper, although the reasons for the metallicity 
tendencies were properly clarified. The main concern about 
finding a metal in the bulk is related with issues of Anderson 
localization, particularly in 2D, as discussed in this appendix. 

In classical transport theory, which relies on weak scatter- 
ing, the conductance G of a cZ-dimensional hypercube of vol- 
ume L d is related with the conductivity a through G=aL d ~ 2 , 
where L is much bigger than the mean free path, I. How- 
ever, Abrahams et. al. have shownS2i - based on renor- 
malization group ideas and perturbation theory - that there is 
no mobility edge in 2D, hence the conductance should van- 
ish as L^oo, where quantum interference leads to Anderson 
localization. According to this picture, the correct scaling of 
the conductivity in one, two, and three dimensions is given 
byM.: a lv (L)=a -*£(L- I), a 2 n(L)=a - ln(f ), and 

V3d{L)=<jq — jr^il ~t)< respectively, where ctq = ne 2 r jm 
is the Drude conductivity at scale I and r = 1/frkp is the re- 
laxation time. In 3D, it seems possible to have a non-zero 
conductance for macroscopic lattice sizes, but not in a lower 
dimension. 

To compare these formulas for localization with our simu- 
lations, we proceed as follows. In 2D and 3D, the conductivity 
attains its maximum value at the disorder strength A=0.8 and 
A=0.9, respectively. We define the difference between con- 
ductivities for two lattices of sizes L\ and L 2 as: 

ACtib^) - ^ - 8 (L2)-a 2 ^°- 8 (L 1 ), (16) 
Ao*g(L 1} L 2 ) = ^r - 9 (i 2 )- O ^=°- 9 (L 1 ), (17) 

where "toy" indicates that these are the numerical results for 
the toy model. The corresponding "localization" formulas are 
given by: 

Aa*g(Li,L2) = -ln(^i), (18) 
Aa^(L 1 ,L2) = -l(-L--L). (19) 

In all of the above, the units for conductivity are "e 2 //i". The 
comparison between these formulas is summarized in Table|l] 
for the 2D case. 

TABLE I: Differences in conductivities between pairs of lattices in 
2D, using data from both the toy model and localization formulas. 



L2 Li AgffifLi.La) A^° D c (Li,L2) 

24 20 -0.1167 -0.1161 

24 16 -0.1844 -0.2581 

24 12 -0.3735 -0.4413 

20 16 -0.0677 -0.1421 

20 12 -0.2568 -0.3252 

16 12 -0.1891 -0.1831 



The similarity of the numbers leads to the conclusion that 
the numerical results for the 2D toy model in the "metallic 
regime" could still lead to an Anderson insulator as the lattice 
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size grows. However, obviously the mechanism discussed in so severe in 3D, giving confidence that the state stabilized at 
this paper is mainly to be applied to 3D manganites (or bilay- intermediate A's is indeed a metal in the bulk limit, 
ered systems). It is known that Anderson localization is not 
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